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We calculated in detail the angular distribution of 7-rays and the resulting shape of the 7-ray line 
produced by the nuclear deexcitation of the 4.439 MeV state of '^^C following proton and a-particle 
interactions with ^'^C and ^^O in the energy range from threshold to 100 MeV per nucleon, making 
use of available experimental data. In the proton energy range from 8.6 to 20 MeV, the extensive 
data set of a recent accelerator experiment on 7-ray line shapes and angular distributions was used to 
deduce parametrizations for the 7-ray emission of the 2+, 4.439 MeV state of ^^C following inelastic 
proton scattering off ^'^C and proton induced spallation of ^^O. At higher proton energies and for 
■ a-particle induced reactions, optical model calculations were the main source to obtain the needed 

' reaction parameters for the calculation of 7-ray line shapes and angular distributions. Line shapes 

are predicted for various interaction scenarios of accelerated protons and a-particles in solar flares. 
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I. INTRODUCTION 



The 7-ray line at 4.438 MeV from the deexcitation of the first excited state of ^^C is one of the strongest nuclear lines, 
J — ' clearly visible in several energetic solar flares, which were observed in particular by the Gamma Ray Spectrometer 
[ — ' aboard the SMM satellite [Q and one of the best candidates for an observation at 7-ray energies of the interaction 
C\j 1: of low energy cosmic rays in nearby molecular clouds Q . It is mainly produced by inelastic scattering of energetic 
■ protons and a-particles off ^^C and by the spallation of ^^O by the same particles This line is also produced in 
1; reverse kinematics by accelerated ^"^C and ^^O bombarding the ambient hydrogen and helium nuclei, with however 
a large Doppler broadening due to the high velocities of the emitting ^^C nucleus, which makes its detection much 
more difficult. In the following we will refer to 7-ray lines produced by light ion bombardment of the ambient gas as 
i-S^ ■ the narrow component and to 7-rays produced by accelerated heavy ions as the broad component. 
Jp-i Accordingly, the intensity and shape of the narrow component of the 4.438 MeV line depend on the properties 
Q of the accelerated light particle spectrum and the ^^C and ^^O abundances in the interaction site. The narrow line 
intensities from various isotopes can be used to derive accelerated particle and ambient gas abundances, as it has been 
^ done for example by Murphy et al. Q for the solar fiare of 1981, April 27. More detailed information on the energy 
spectrum of the accelerated particles and their directional distribution at the interaction site may however only be 
obtained by a detailed line shape analysis. 

Murphy, Kozlovsky and Ramaty ^ investigated the impact of different accelerated particle angular distributions 
^ , on the narrow component of the 4.438 MeV 7-ray line produced by the p + ^^C reaction. The shape of the broad 
component has been discussed by Bykov, Bozhokin and Bloemen Q and by Kozlovsky, Ramaty and Lingenfelter 
0. These studies are based on the experimentally well-known excitation functions, some available experimental data 
on differential cross sections, few published results on measured line shapes, and theoretical arguments. For the 
spallation reaction, the situation is similar with respect to available experimental data. The excitation function is 
experimentally determined from threshold up to over 100 MeV, but few data are available on differential cross sections 
and line shapes. 

In order to put the calculations for the 4.438 MeV line on a firmer basis, we use additionally to the already published 
data the extensive data set on 7-ray angular distributions and 7-ray line shapes for the ^^C(p,p7) and -^^0(p,pa7) 
reaction, obtained in a recent accelerator experiment at Orsay [||. For the spallation reaction, the reproduction of the 
observed line shapes for proton energies between 16 and 20 MeV proved to be only sensitive to the mean excitation 
energy of the a-particle emitting intermediate states of ^^O. Estimations of the mean excitation energy for higher 
proton energies and for the differential cross sections are based on published data in the proton energy range 40-50 
MeV and optical model (OM) calculations. 

For proton inelastic scattering, we adopted a method similar to that of Murphy, Kozlovsky and Ramaty for line 
shape calculations. Two independent amplitudes for the magnetic substate population of the 2+ level and a phase shift 
were adjusted to reproduce both the experimental line shapes and the laboratory 7-ray angular distributions in the 
proton energy range from 8.6 MeV to 19.75 MeV. The extrapolations to higher energies are based on OM calculations 
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using phenomenological potential parameters, which proved to give an excellent reproduction of measured differential 
cross sections for nucleon elastic and inelastic scattering off ^^C at projectile energies E « 20 - 100 MeV ||]. 

a-particle inelastic scattering off ^^C and spallation of ^^O suffer from a complete absence of published data on 
line shapes. Therefore, estimations of the mean excitation energy in -'^^O for the spallation reaction were entirely 
based on OM calculations and an experimental a-particle spectrum at Ea — 60 MeV. For the inelastic scattering 
reaction, differential cross sections da/dfla and some laboratory 7-ray angular distributions have been used to find 
the parameters for the 7-ray emission based on the population amplitude method at a-particle energies below 20 
MeV, while at higher energies the needed parameters were obtained from OM calculations. 

In the last section we discuss the production of the 4.438 MeV line in solar flares and show some calculated shapes 
for the narrow component of the 4.438 MeV line induced by proton and a-particle interactions in the solar atmosphere 
as it may be observed by the HESSI spacecraft pO[ |. 



II. THE SPALLATION REACTION ^^0{p,x-/y^C 

Total cross section data are reported by Dyer et al. for Kj ~ 14-23 MeV, Lesko et al. for Ep — 20, 30, 33, 
40 and 50 MeV, Lang et al. ^ a.t Ep = 40 MeV, Zobel et al. Q at £'p = 12.1, 28.2, 48.3 and 145 MeV and by Foley 
et al. jl5| at Ej, ^ 146 MeV. From threshold to 23 MeV, we used the data of Ref. [|ll|, which are in agreement with 
those of Refs. |12| and |0. There is, however, some disagreement between the different data in the energy interval 
Ep = 40-50 MeV. The cross section values at Ep = 40 MeV of Ref. ||l^ and Ref. Q differ by more than a factor of 
two. The same holds for the values around 50 MeV of Ref. |Q and Ref. |0. As no systematic deviation of one of 
the data sets can be observed, no cross section value was discarded. We therefore fitted the excitation function curve 
above 23 MeV. Assigning the same weight to each data point, the overall shape of the excitation function curve from 
23 MeV to 146 MeV could be best fitted by a power law plus a constant. The fit is very close to the curve in Fig. 3 
of Ref. 1^ for this reaction. The obtained cross sections values are given in Table |. 

For proton energies below 50 MeV, that reaction should proceed principally from a sequential process with the 
excitation of intermediate states of ^^O, decaying subsequently to a -I- ^^€4 439. Reactions with other particles than 
the a-particle in the outgoing channel, like '^He -I- n or '^H -I- p have thresholds « 20 MeV higher, around Ep — 35 
MeV, and should be negligible for the production of ^^€4 439 in solar flares. The sequential character can be deduced 
from a study of the ^^0(p,pa)^^C reaction at a proton energy of 46.8 MeV, where no evidence for quasi-elastic p-a 
scattering was found ||l^. Still other reaction paths, such as ^®0(p,a)^'^N* ^ p + ^^€4 439 were also found negligible. 
At lower energies these reaction mechanisms should be still weaker. We adopted in the following exclusively the 
sequential process ^^0(p,p')^^0* ^ a + ^^04 439. 

Parameters entering into the line shape determination are therefore the relative contribution of the different excited 
states in ^^O and for each excited state the differential cross section with respect to the recoiling ^^O* scattering 
angle, the ^^C*-"'^^0 angular correlation and the angular correlation. Given the scarcity of experimental data 

and the number of excited states involved, it is highly impracticable if not impossible to undergo a full and correct 
calculation for numerous proton energies from threshold to 100 MeV. One can, however, reasonably assume that the 
average angular distribution of the a-^^C* decay in the recoiling ^^O* system and the angular distribution of the 
7-ray emission in the ^^C* system are isotropic, because of the contribution of many ^^O levels with different spin 
and parity (see the discussion in Ref. [p^ ). 



A. Ep = 16 - 20 MeV 



In this energy range, experimental line shapes are available from a recent accelerator experiment at Orsay [g[ at 
six different 7-ray detection angles, ranging from 45° to 145° with respect to the incoming beam direction and at 
four proton energies. These experimental line shapes were compared with calculated ones obtained with the use of a 
Monte-Carlo type program. In this program reaction parameter inputs were the differential cross section with respect 
to the proton center-of-mass scattering angle and the mean excitation energy in ^^O. As mentionned above, the vector 
of the a-^^C* decay in the ^^O* system and the 7-ray emission in the ^^C* system were taken isotropically. Slowing 
down of the carbon in the target with exponentially decaying 7-ray emission probability corresponding to the lifetime 
of the 4.439 MeV state, as well as the full detector geometry and resolution were taken into account to generate 
theoretical detector spectra. Stopping powers were taken from the code TRIM pTf . 

Optical-model calculations of inelastic proton scattering off ^^O to states from 12.75 MeV to 17 MeV, with optical 
potential parameters from the compilation of Percy and Percy [ p^ were performed to obtain the differential inelastic 
scattering cross sections da/dQp^^. The calculated cross sections were found to be roughly isotropic in this proton 
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energy range. As furthermore the calculated line shapes proved to be not very sensitive to the details of the differential 
cross section, we simply used a constant dcr/dfip^^ . The mean excitation energy in ^^O was then the only parameter 
that has been adjusted to reproduce the experimental 7-ray spectra. 

Figure |l| shows the result of a line shape calculation at Ep — 20 McV for the six detection angles of the experiment. 
Position, line shape and relative intensities of the experimental lines are fairly well reproduced by the calculation, 
which gives credit to the validity of the assumption of essentially isotropic angular distributions. We also checked that 
calculations with a Gaussian distribution of the excitation energy centered at £x resulted in practically identical line 
shapes. Those calculations have been made at the four proton energies where experimental spectra were extracted. 
The obtained mean excitation energies are presented in Table ffl. 



B. Ep = 20-100 MeV 



In this proton energy range, some experimental information on line shapes and on inelastic cross sections of excited 
levels in ^^O are available at Ep = 40, 45 and 46.8 MeV [|6). Lang et al. ^ deduced relative contributions 

to the ^^€4 production for ^^O states ranging from E^ — 12.53 MeV to 28 MeV from the cross sections at Qp = 
40° and Ep = 45 MeV of Buenerd et al. ||l8|. For excitation energies above 28 MeV, virtually no contribution has 
been observed in the i60(p,py'^0 reaction at Ep = 45 McV and in the ^^0{p,pay^Cl 439 reaction at Ep = 46.8 
MeV 

We used as the basis of our calculations eight representative states Ei for the ^^O levels, with the relative contribution 
to the ^^04 439 production Pf. 

^ £k^(40°) (1) 



given by Lang et al. |13| in their table I for Ep = 45 MeV. Here, F^^ and FJ^j mean the partial decay width to the 
a + ^^CX 439 channel and the total width of the state i, respectively. In detail, we chose the following representative 
states with excitation energy Ei: 

(1) El = 12.75 MeV, 2+, representing the 12.53, 13.02 and 13.26 MeV levels of Ref. Jp). 

(2) E2 = 15.1 MeV, 2+, representing the 13.97, 15.27 and 15.50 MeV levels of Ref. |3|| 

(3) i?3_8 = 17, 19, 21, 23, 25 and 27 MeV, 2+, representing the respective excitation energy range E^ — Ei ± 1 
MeV of Ref. (l|. 

The choice of exclusively 1=2 excitations to 2+ states has two reasons. First, calculations of inelastic proton 
scattering to the giant resonance region at E^ = 21-28 MeV indicate good agreement of 1=2 excitation with the 
experimental data |Q . Second, the dependence of the cross section on the proton energy proved to be not very sensitive 
to the multipolarity, and the angular dependence of the cross section for 1=2 excitations is somehow intermediate 
between 1=1 and Z=3 excitations, simulating fairly well an averaged angular distribution for the three multipolarities. 
Non-natural parity states and excitations with Z>3 are believed not to change much the overall trends. 

We also assumed that the total cross section for each state is proportional to the differential cross section at 40°. 
This was supported by OM calculations for Ep = 45 MeV, where the ratio of the inelastic scattering cross section ai 
to the differential cross section dai/dil at 40° was found constant to within 20% for the different excited states. 

The mean excitation energy for Ep = 45 MeV is then simply defined as: 

= (2) 

To obtain the mean excitation energy at other proton energies, we made extensive OM calculations with the code 
Ecis94 20 1, using the optical potential parameters of the compilation of Percy and Percy [M. For each proton energy 
we calculated the inelastic scattering cross section for the eight representative states of O. All calculations were 
restricted to direct one-step excitation of the respective state in -'^^O, which should be the main excitation mechanism 
for the states of interest at these proton energies. 

Calculations were done for the following proton energies: Ep = 22.5, 25, 30, 40, 45, 66 and 100 MeV. The mean 
excitation energy as a function of proton energy was calculated in the following way, 

o (r,^ P^ a^{Ep)/a.{45 MeV) 
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The differential cross sections have been obtained in the same way by averaging over the eight representative states. 
In order to facihtate hne shape calculations for a wide range of proton energies, we parametrized the differential cross 
section by a simple function: 

ailp 

where O is the proton diffusion angle in the center-of-mass system. It was checked that this simple function produced 
essentially the same line shapes as the differential cross sections from the OM calculations. The values of £x and p 
are presented in Table |. 



III. THE INELASTIC SCATTERING REACTION ^^C{p,p'y) 

The cross section excitation function ^{Ep) was taken from the measurements of Dyer et al. |lT| from threshold 
to 23 MeV, above that from the data of Lang et al. and Lesko et al. up to 85 MeV, then from the curve of 
Ramaty, Kozlovsky and Lingenfelter |^ up to 100 MeV. Cross section values are reported in Table ||. 

Production of 4.438 MeV 7-ray emission by inelastic scattering off ^^C originates mainly from direct excitation of 
the 2+, 4.439 MeV level. All higher- lying levels of ^^C are particle-unbound and have small to very small 7-branching 
ratios to the 4.439 MeV state |2|], making 7-ray cascade contributions negligible. The 7-ray line shape and the 
7-ray angular distribution are in this case completely determined by the double differential inelastic cross section 

a / dflpdflj . It can be expressed as a product of the differential cross section da/dflp and the angular correlation 
function W (see for example Satchler ||2^): 

- ^i:WiVk'rk,) (5) 



dflpdQy dflp An 

where k and k' are the wave vectors of the incoming and outgoing proton, respectively, and fc-y is the vector of the 
emitted 7-ray. 

Although a wealth of differential cross sections are available for this reaction, proton-7 correlation measurements 
exist only for few proton energies. The angular correlation function must therefore be extracted either from experi- 
mental data, for example by fitting of 7-ray angular distributions and line shapes or be obtained from OM calculations. 
The excitation function for the 4.438 MeV 7-production cross section exhibits pronounced fluctuations from threshold 
up to w 15 MeV (see Ref. |^ or Ref. |l^). This indicates strong contributions of compound nucleus resonances ^^C-|-p 
i3p^* p_|_12q* ^Yie inelastic scattering. 

Because OM calculations are only suited for the direct reaction mechanism, their use in this energy range was 
discarded. We used the measured line shapes and 7-ray angular distributions of the Orsay experiment |^ and published 
data for da/dilp for Ep = 8.6 - 20 MeV to adjust the parameters of the correlation function. Some available 7-ray 
angular distributions of the Washington experiment [ll| were also used for comparison. The parameter adjustment was 
achieved by reproducing as closely as possible the experimental data in extensive calculations with the Monte-Carlo 
type program simulating the reaction and the detection setup of the Orsay experiment. 

At higher proton energies, Meigooni et al. |^ provide a phenomenological optical potential for neutron and proton 
scattering off ^^C, which reproduces fairly well elastic and inelastic scattering angular distributions. The OM calcu- 
lations were used to obtain the necessary input parameters for the calculation of the double differential cross section 
dV/dOpdr?^. 



A. Ep = 8.6 - 20 MeV 

We adopted the formalism used in Refs. |^ and ||^ for the proton-7 angular correlation and the calculation of 
the 7-ray line shape. In this formalism, the angular correlation function W is expressed in the rest frame of the 
recoiling excited ^^04 439 nucleus, whose recoil angle and velocity is kinematically fixed by the incoming and outgoing 
proton wave vectors k and k'. The correlation function W then reduces to the 7-ray angular distribution W{&, $) in 
the carbon rest frame. It is proportional to the radiation pattern dP/dVL, which depends on the magnetic substate 
population amplitudes am{k, k') of the 4.439 MeV, 2+ state of ^^C (see e.g. Jackson p3|): 

^ = ^ I E *«™(fc,fc')^2™(e,$)p (6) 

m=-2 
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where the are the vector spherical harmonics for quadrupole radiation, and 0, $ the polar and azimuthal angles 
of the 7-ray emission in the rest frame of the excited carbon. The aim was to find a set of population amplitudes 
am{k,k') for each proton energy of the Orsay experiment, which gives simultaneously a good reproduction of the 
experimental line shapes and the laboratory 7-ray angular distribution. 

Kozlovsky, Ramaty and Lingenfelter obtained reasonable fits of the six experimental line shapes at Ep — 23 
MeV from ref. by adjusting only two independent real amplitudes, which furthermore were independent of the 
proton diffusion angle. We started our search with the same simple basis and adjusted the two independent relative 
amplitudes a2/aa and ai/a^ (with a„j real, a_m = Om and X^'^m ~ l)i independent of the carbon recoil angle. 
However, these two parameters were not sufficient and it proved to be necessary for a simultaneous reproduction of 
the line shapes and the 7-ray angular distribution to introduce a third independent parameter. We found that a phase 
shift A<i> for the azimuthal angle improved considerably the fits and gave a satisfactory reproduction of the Orsay 
data at most proton energies. 

Practically, for the search of the population amplitudes, we used as in the case of the spallation reaction a Monte- 
Carlo type program simulating the nuclear reaction and the detection system of the Orsay experiment. For each proton 
center-of-mass diffusion angle, the ^^C\ 439 recoil direction and energy were calculated in the laboratory system. For 
the 7-ray emission, we chose a coordinate system with origin in the recoiling carbon nucleus, the z-axis perpendicular 
to the scattering plane and the x-axis in the direction of the carbon velocity vector. O defines the angle between the 
7-ray vector and the z-axis, while $ is defined as the angle between the x-axis and the projection of the 7-ray vector 
on the scattering plane, counted counterclockwise from the x-axis. Subsequently, the 7-ray vector was calculated in 
the laboratory system and the energy of the 7-ray was stored when it was in the solid angle of one of the detectors. 
Again, slowing down of the recoiling excited -'^^C in the target and decay of the 4.439 MeV state with half- life T1/2 = 
42 fs were taken into account. 

Because of the presence of multiple compound nucleus resonances in this energy range, calculations were done 
in narrow proton energy steps between 8.6 and 20 MeV. For each energy, population amplitudes were adjusted to 
reproduce the 7-ray line shapes of the six detectors and special care was taken to obtain a simultaneously good 
reproduction of the laboratory 7-ray angular distribution. For da/cWlp, experimental cross sections of Peellc |^ for 
£'p = 14 - 19.4 MeV and of Barnard, Swint and Clegg ^ for Ep = 8.6 - 11.6 MeV were used. For Ep = 11.6 - 14 
MeV we interpolated between the differential cross sections of Ref. at 11.6 MeV and Ref. at 14 MeV. 

The results of the adjustments are presented in Table || and an example of calculated line shapes and laboratory 
7-ray angular distribution is shown in Fig. p| together with data of the Orsay H and Washington ||ll| experiment. 



B. Ep = 20 - 100 MeV 

Experimental data in this energy domain are relatively scarce, and no systematic measurements of differential cross 
sections da/dVlp and da/AQ,^ or line shapes are available. However, Meigooni et al. j^] performed a systematic study 
of nucleon elastic and inelastic scattering off ^^C and found an energy dependent optical potential parameter set, 
which provides good fits to measured differential cross sections in a wide angular range and for energies above the 
pronounced compound nucleus contributions (« 15 MeV) up to approximately 100 MeV. We therefore based our 
calculations entirely on OM calculations with their parameter set. 

These calculations were done with the program Ecis94 to obtain the scattering amplitudes TM'^,M'^,MAMa i^, k') 
necessary for the construction of the correlation function W (Eq. ^). Ma {M'j^) and Ma {M'^) are the z-components of 
the ^^C and the proton spin in the incoming (outgoing) channel, respectively. Choosing the z-axis along the incoming 
proton direction k, the amplitudes, correlation function and consequently the proton and 7-emission angle are directly 
given in a space-fixed system with respect to the beam direction of the experiment. The correlation function W for 
emission of the 7-ray in the direction k^ is given by: 

W{k,k',k^) = tkq{k,k') Rig (7) 

kq 

where t^q denotes the polarization tensor. It is constructed from the scattering amplitudes T and contains all 
information on the state of polarization of the excited carbon after scattering (see Satchler section 10.3.3). Rkq 
are the radiation tensors: 

Rkq = Rkh) \I^^Y,\e^,<f^) (8) 

where Rk{l) are the gamma radiation parameters and Y^f the spherical harmonics. 
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Introducing this correlation function in a program similar to the Monte-Carlo type program described above, 
we calculated theoretical line shapes and 7-ray angular distributions in the laboratory, using the differential cross 
sections da/dSlp obtained from OM calculations. Comparison with the Orsay data for Ep = 16-20 MeV gave fairly 
good reproductions of the line shapes and 7-ray angular distributions. An example is shown in Fig. ^ for Ep — 19.75 
MeV. Given the tendency that the differential cross sections are better reproduced with increasing proton energy (see 
Ref. [^), the calculated line shapes and 7-ray angular distributions should be fairly realistic above 20 MeV. 

Those calculations are, however, quite complex, involving for each proton energy the input of five complex scattering 
amplitudes depending furthermore on the proton center-of-mass scattering angle. In order to facilitate line shape 
calculations for proton spectra covering a wide energy range as in solar flares or cosmic rays, the line shapes and 
7-ray angular distributions from the OM calculations for Ep > 20 MeV, were adjusted with the population amplitude 
method described in subsection III. A. The population amplitudes and phase shifts from these adjustments are reported 
in Table ||. 

IV. ALPHA-PARTICLE INDUCED REACTIONS 

As in the case of proton induced reactions, the main source of 4.438 MeV 7-rays by accelerated a-particles is 
inelastic scattering off ^^C and spallation of ^^O. The cross section excitation functions being quite similar for both 
projectiles, reactions with a-particles should generally not contribute significantly to the 7-ray line due to their low 
abundance, e.g. an observed a/p ratio of 0.035 and 0.0076 in the energetic particle spectrum of two solar flares pTj . 

However, a/p ratios in excess of 0.1 for the accelerated particle spectrum were shown to be necessary for an 
acceptable fit of the 7-ray spectrum of the 1981 April 27 flare a-particle reactions are furthermore favored in the 
case of a very soft particle spectrum, because the thresholds for inelastic scattering off ^^C and spallation of ^^O are 
lower for a-particles (~ 2 and 5 MeV per nucleon, respectively) than for protons (« 6 and 14 MeV, respectively). 
In such cases, an important fraction of the 7-ray production may be due to reactions with a-particles. Therefore, a 
treatment of a-particle reactions similar to the above described studies on proton reactions seems worthwhile. 

Unfortunately, experimental line shapes are practically absent for a-particle inelastic scattering off ^^C and do not 
exist for a-particle induced spallation of ^^O. Concerning spallation of ^^O, this may not be very problematic because 
of the excitation of many intermediate levels in ^^O resulting in essentially isotropic 7-ray angular distributions as for 
^^0(p,pa7). Optical model calculations should then provide as for proton induced spallation a reasonable estimation 
of the mean excitation energy in ^^O and approximate differential inelastic cross sections. On the contrary, inelastic 
a-particle scattering off ^^C is certainly dominated at energies below Ea ~ 30 MeV by compound nucleus resonances 
due to the pronounced a-particle structure of ^^C and ^^O, making OM calculations very hazardous in this energy 
range. 

A. The spallation reaction ^®0(q,x7)^^C 

We made the same assumptions as for the case of proton induced spallation of ^^O. Only the sequential process 
of excitation of intermediate levels in ^^O with subsequent decay in a and ^^04439 was considered. The angular 
distributions of the excited carbon in the system of the excited oxygen and of the 7-ray in the system of the excited 
carbon were both considered isotropic. 

Because of the absence of experimental line shapes, the estimation of the mean excitation energy in ^^O was entirely 
based on OM calculations and the a-particle energy spectrum of Lang et al. |jl^ , obtained at Ea — 60 MeV and at 
a scattering angle of 12.5°. Contrary to proton scattering, the ratio of the total inelastic scattering cross section to 
the differential cross section at 12.5° shows a systematic dependance on the excitation energy in the considered range 
Ex = 12.53 - 28 MeV. Therefore, the relative contributions dai/dQ{12.5°) extracted from the spectrum of Ref. p3[ , 
were multiplied by this ratio to obtain the relative inelastic scattering cross sections ct; at Ea = 60 MeV. Finally, for 
the partial 7-production probability at 60 MeV, the inelastic scattering cross sections were multiplied with the 
branching ratios F^yrj^j, proposed by Ref. 

Extensive OM calculations were done for a-particle energies Ea = 18 - 400 MeV and the same eight representative 
excited "'^^0-states as in subsection II. B to estimate their relative contribution at other projectile energies. As in the 
proton case, we restricted the calculations to direct one-step I — 2 excitations. The energy-dependent phenomenolog- 
ical optical potential of Michel et al. |2^ , which gives good to very good fits of elastic differential cross sections from 
about 20 MeV to 146 MeV, were used in the entire energy range. From the same calculations, averaged differential 
cross sections were constructed and their angular distributions fitted by a parameterization as in subsection II. B, but 
with the addition of a term cos^(O) to take account of the backward enhancement of the differential cross sections: 
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^ = Cie-f^ + C2Cos4(e) (9) 

For this reaction, the 7-ray production cross section has only been measured at five a-particle energies by Dyer et 
al. H] from 22 to 26 MeV and by Zobel et al. [|l|| at i;^, = 52 MeV. An estimated cross section curve is given in Ref. 

, however based on the then only available cross section data of Ref. and limited to a-particle energies below 
20 MeV per nucleon. We therefore made a new estimation extending the energy range from threshold to 100 MeV 
per nucleon. Our best guess, based on the energy dependence of the (a,4pxn) spallation data of Lange et al. and 
the OM calculations is a linear interpolation between the cross sections at 26 MeV and 52 MeV. Below 22 MeV and 
above 52 MeV, we take the energy dependence of the OM calculations. 

The results of the calculations for the mean excitation energy, the inelastic scattering angular distributions and the 



7-ray production cross sections are presented in Table [I] 



B. The inelastic scattering reaction ^^C{a,aj) 

As already mentionned, practically no line shapes are available for this reaction and OM calculations at low a- 
particle energies are not believed to be very realistic. To get anyway an idea of the magnetic substate population 
amplitudes, we employed two different methods at energies below 30 MeV. From Ea = 8-17 MeV, we used the 
systematic measurements of differential cross sections da/d^a and laboratory 7-ray angular distributions da/dQj 
of Mitchell, Carter and Davis and Ophel et al. We then used the Monte-Carlo type program described in 
subsection III. A to find a set of population amplitudes am (with am—a^m, and here throughout a±i=0 for 0+ 2+ 
excitations induced by inelastic a-particle scattering on spin-0 targets like ^^C) and phase shifts A$ which reproduces 
the measured 7-ray angular distributions. 

Morgan and Hobbie measured in small energy steps between Ea = 19-30 MeV inelastic differential cross sections. 
However, no 7-ray angular distributions or line shape measurements are available in that energy range. Therefore, 
we used OM calculations despite the fact that compound nucleus resonances are probably important below E^ — 30 
MeV. In order to reduce the effect of individual compound nucleus resonances, we averaged several differential cross 



sections of Ref. |33| centered around the a-particle energies of 21 and 25 MeV. Between 30 MeV and 200 MeV, we 
used the inelastic scattering data of Burdzik and Heymann at 32.5 MeV of Baron, Leonard and Stewart at 
41 MeV, of D'Agostino et al. at 90 MeV of Smith et al. at 139 MeV and of Tatischeff and Brissaud at 166 
MeV l38|- With these differential cross sections OM fits were then done, again with the optical potential of Michel 
et al p8| . Results of the OM fits in the energy range Ea — V! - 166 MeV are shown in Fig. ^. Above 200 MeV, no 
inelastic differential cross sections are published. We continued to use the optical potential shape of Michel et al. for 
the inelastic scattering calculations, extrapolating up to 400 MeV the optical potential found at 166 MeV by following 
the energy dependence of the potential parameters given by Michel et al. for a-particle scattering off ^^O. 

Again, as for proton inelastic scattering off ^^C, the line shapes and 7-ray angular distributions predicted by the 
OM calculations were used to deduce population amplitudes and phase shifts for a consistent parameterization of the 
7-ray emission in proton and a-particle inelastic scattering off ^^C in the energy range of interest for solar flares and 
low energy cosmic ray s. The cross section excitation function is taken from Dyer et al. |2^] for Ea = 8-26 MeV, from 
the data of Refs. [p4[-|3q] at 32.5, 41, 90, 139 and 166 MeV, respectively, and from our OM calculations at 250 and 



400 MeV. Results for the population amplitudes, phase shifts and total cross sections are presented in Table IV 



V. PREDICTED LINE SHAPES IN SOLAR FLARES AND DISCUSSION 

With the above deduced reaction parameters, the line shape of the 4.438 MeV 7-ray emitted from the first excited 
state of ^^C for proton and a-particle interactions with ^^C and ^^O can be calculated for energies from threshold up 
to 100 MeV per nucleon. To illustrate the line shapes that can be expected from solar flares, we use a simple model of 
7-ray production by energetic particle interactions in a thick target. For the accelerated particle energy distribution 
N{E), we choose a spectrum as it results from diffuse shock acceleration [p7| : 

N{E) oc - p-" exp{-E/Eo) (10) 

V 

where v and p are the particle velocity and the momentum per nucleon, respectively; s and Eq are the s pec tral 
index and cutoff energy, specific to a given acceleration site. The energies labelled E, Eq, Emax, E' in Eqs. (10|Tl|) 



and in the following are expressed in energy per nucleon. The total number of 4.438 MeV 7-rays produced by 
interactions of projectiles of type i (p or a) with target nuclei of type j is then given by: 
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where [Aj] is the abundance of isotope j with respect to hydrogen and n the density of hydrogen atoms at the 
interaction site; aij{E) is the 7-ray production cross section and dE/dxi{E) is the stopping power of particle i. We 
made the simphfying assumption that the energy loss is only due to the electronic stopping power in a neutral gas, 
consisting of hydrogen and helium with [He\/[H] = 0.1. The Bethe formula was used for the energy loss calculation. 

We adopted two sets of s and E^: s = 3.3; Eq = 30 MeV and s ~ 2.4; Eo > 300 MeV corresponding to the fit of 
interplanetary proton observations of two distinct solar flares We set [O] — 0.068%, [C] = 0.042% as in Ref. [ p7| , 
Emax = 2 X Eq and an a/p-ratio of 0.1. 

In both cases, proton interactions dominate; they account for 82% and 84% of the 7-ray production with the softer 
and harder spectrum, respectively. With the softer spectrum, the 7-ray production is dominated by proton inelastic 
scattering off ^^C (57%), of which in particular 85% of the 7-rays are produced by interactions at proton energies 
below 20 MeV. In order to get an overall impression of the quality of reproduction of the 7-ray lines from the important 
low energy proton inelastic scattering reactions, energy integrated data from the Orsay experiment were produced 
and compared with the calculations. Line shapes and 7-ray angular distributions for energy integrated experimental 
and calculated data between 8 MeV and 20 MeV for proton inelastic scattering off ^^C, are represented in Fig. ||. 

For the relatively hard spectrum with s = 2.4 and setting Eq = 300 MeV, proton induced spallation of -'^^O is the 
strongest channel, accounting for 46% of the 4.438 MeV 7-rays. Here, proton interactions at energies below 20 MeV 
are responsible for 24% of the total 7-ray production, a and proton interactions at energies below 100 MeV account 
for 83% of the produced 7-rays. These values, however, have to be regarded as lower limits only, because secondary 
protons of low and intermediate energy produced by high energy proton and a-particles and escape of accelerated 
particles into interplanetary space are not taken into account. 

For the prediction of 7-ray line profiles from solar fiares we calculate the interaction probability as in Eq. (^ij). No 
absorption or Compton scattering of the produced 4.438 MeV 7-ray in the solar atmosphere is taken into account. 
Three different energetic particle angular distributions have been investigated, similar to the distributions investigated 
by Murphy, Kozlovsky and Ramaty |^]: an isotropic angular distribution, a fan-beam fiare and a downward flare. 
The narrow component of the 4.438 MeV 7-ray is shown in Figs. ^, ^and|| for the different energetic particle angular 
distributions and energy spectra, all with an a/p-ratio of 0.1. All curves are arbitrarily normalized to a total kinetic 
energy in proton and a-particles with E >1 MeV of 10* erg. 

The narrow component is clearly dominated by the p-|-^^C reaction for the soft particle spectrum. Comparing curves 
1 and 2 of Fig. ^, which shows the predicted line shapes for isotropic incident particle spectra, it seems possible to 
estimate the spectral hardness by the width of the line and the characteristic drop around 4.438 MeV only visible 
for the hard spectrum. This drop, also seen in curve 5 of Fig. is not predicted by other parameterizations, as for 
example in Ref. j|] . 

Figure shows the line profiles resulting from a downward directed energetic particle distribution at the solar 
limb. The parameterization of Ref. |^ for the population amplitudes reproduces qualitatively the characteristic line 
shape with two intensity maxima above and below the minimum at 4.438 MeV, but it underestimates slightly the 
maximum-to-minimum ratio and overestimates the intensity (by wl0%) emitted at 90°. Line-shapes produced by 
fan-beam distributions and an isotropic distribution as in Fig. g| will be difficult to distinguish, showing only small 
differences in the overall line profile. 

In summary, with the extracted nuclear reaction parameters in this work, line shape calculations for proton and 
a-particle induced reactions on ^^C and ^^O can be done for all possible energetic particle angular distributions and 
for particle energies from threshold to about 100 MeV per nucleon, which is the important energy domain for solar 
fiares and low energy cosmic rays. These parameters can also be straightforwardly used for the inverse reactions, 
energetic ^^C and ^^O bombarding proton and a-particles. 

It is the first parameterization for line shape calculations, which is largely based on experimental data and optical 
model calculations. Especially, the calculation of 7-ray line shapes from proton induced reactions is certainly very close 
to reality, as illustrated by the good reproductions of experimental line shapes in Figs. |^ and ^; at higher energies, the 
optical potential of Ref. |^ is believed to provide line shape reproductions of the same quality. Only line shapes from 
a-particle induced inelastic scattering below w 10 MeV per nucleon have to be regarded with some caution because of 
the only moderately successfull OM fits of differential inelastic cross sections. Here, new experimental data for 7-ray 
line shapes and angular distributions could greatly improve the parameterization. For a-particle induced spallation 
of ^^O above Ea = 26 MeV, 7-ray production cross sections are needed. 

This study confirms qualitatively older parameterizations for proton induced reactions, constructed with much less 
experimental data input but some of its predictions of line-profile details are slightly different. It shows, that 

extreme particle angular distributions, like isotropic or unidirectional ones can be nicely characterized, with even 
some hope to get an idea of the spectral hardness of the accelerated particle spectrum. High resolution detectors. 
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like the Ge detectors on the HESSI spacecraft jl^] and on the INTEGRAL observatory could be able to resolve 
the fine structures in the line profiles which carry this information. Together with other 7-ray lines, one may then be 
able to extract most flare parameters, like the spectrum, composition and the angular distribution of the accelerated 
particles, as well as the isotopic composition of the ambient medium. 
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TABLE I. Cross section a and mean excitation energies for the spallation reaction ^^0(p,pa)'^^C4.439 and parameter p of 
the inelastic scattering angular distributions ^®0(p,p'). 

{Ep = 16-20 McV): Mean excitation energies £x in ^*'0 from the adjustment of the 4.438 MeV 7-ray lines of the Orsay 
experiment [8]. Isotropic inelastic scattering angular distributions were used (p=0). 

(Ep = 22.5 - 100 MeV): Mean excitation energies £x in and fit-parameter p of the inelastic scattering angular distributions 
^•^©(pip') obtained by OM calculations. For _E'p<16 MeV, the values of 16 MeV were used. 



Ep (MeV) 


a (mb) 


Ex (MeV) 


P (deg ^) 


14 


3 


12.75 


.0 


15 


16 


12.75 


.0 


16 


40 


12.75 


.0 


17 


64 


13.25 


.0 


18 


85 


14.0 


.0 


20 


140 


14.75 


.0 


22.5 


156 


14.7 


.001 


25 


137 


15.5 


.008 


30 


95 


17.1 


.009 


40 


55 


18.0 


.027 


45 


44 


18.3 


.030 


66 


24 


18.5 


.057 


100 


14 


18.5 


.078 
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TABLE II. Cross section a and population amplitudes for the inelastic scattering reaction ^^C(p,p)'^^C4.43g. 
{Ep = 8.6 - 19.75 MeV): Result of the adjustment of the population amplitudes am and the phase shift A$ to measured line 
shapes from the Orsay experiment [8] and 7-ray angular distributions from the Orsay experiment and from the Washington 
experiment [11]. 

{Ep = 25 - 100 MeV): Result of the adjustment of the population amplitudes am and the phase shift A3> to line shapes and 7-ray 
angular distributions resulting from OM calculations. Population amplitudes below 8.6 MeV are an average of the extracted 
values from 8.6-10.6 MeV. 

The amplitudes are normalized such that = 1 and a-m = am- 



Ep (MeV) 


(7 (mb) 


do 


ai 


a2 


A$ (des) 


5 


51 


.64 


.27 


.47 


30 


6 


55 


.64 


.27 


.47 


30 


7 


108 


.64 


.27 


.47 


30 


8 


240 


.64 


.27 


.47 


30 


8.6 


270 


.57 


.28 


.51 


35 


9.0 


290 


.66 


.26 


.46 


30 


9.6 


275 


.71 


.35 


.35 


30 


10.0 


265 


.70 


.28 


.42 


30 


10.6 


296 


.59 


.21 


.53 


20 


11.0 


317 


.55 


.33 


.49 


35 


11.4 


298 


.56 


.34 


.45 


30 


12.0 


270 


.55 


.22 


.55 


25 


12.6 


277 


.42 


.08 


.63 


25 


13.0 


282 


.40 


.12 


.64 


25 


13.6 


281 


.42 


.13 


.63 


25 


14.0 


281 


.42 


.13 


.63 


25 


14.4 


255 


.42 


.13 


.63 


25 


15.2 


217 


.42 


.13 


.63 


25 


16.25 


188 


.33 


.03 


.66 


20 


17.25 


166 


.27 


.03 


.68 


20 


18.25 


153 


.27 


.16 


.66 


15 


19.75 


130 


.29 


.20 


.64 


15 


25 


99 


.25 


.025 


.68 


112.5 


35 


69 


.37 


.0 


.66 


100 


49.5 


23 


.43 


.0 


.64 


100 


75 


14 


.51 


.0 


.61 


95 


100 


11 


.58 


.0 


.58 


90 
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TABLE III. Estimated cross section a of the spallation reaction, fit-parameters p and C2/C1 of the inelastic scattering 
angular distributions ^^0{a,a') and mean excitation energy 8^ of deduced from OM calculations for the spallation reaction 

l«0(a,2Q)12Cl439- 



Ec, (MeV) 


(mb) 


P (deg ^) 




So: (MeV) 


18 


0.22 


.020 


.12 


12.8 


20 


19 


.017 


.01 


12.8 


22 


44 


.012 


-.08 


13.0 


24 


102 


.009 


.06 


13.8 


26 


125 


.008 


.27 


14.2 


28 


124 


.005 


.51 


14.5 


30 


123 


.003 


.61 


14.8 


35 


120 


.001 


1.91 


15.7 


40 


117 


.010 


.94 


16.3 


50 


111 


.019 


.99 


16.6 


60 


139 


.015 


.42 


17.1 


100 


223 


.074 


.0 


18.4 


150 


172 


.125 


.0 


18.8 


250 


100 


.31 


.0 


19.0 


400 


5.'! 


'A 


.0 


19.0 



TABLE IV. Cross section a and results of the adjustment of the population amplitudes am and the phase shift A# to 

measured and calculated line shapes and 7-ray angular distributions for inelastic a-particle scattering to the 4.439 MeV, 2"*" 
state of ^■^C. The amplitudes are normalized such that ain = 1 and a_m = dm,', a±i=0. 



(MeV) 


a (mb) 


Oo 


ai 


a2 


A$ (d( 


8 


48 


.0 


.0 


.71 


-5 


9 


72 


.0 


.0 


.71 


-5 


10 


231 


.0 


.0 


.71 


-5 


11 


347 


.0 


.0 


.71 


-5 


12 


393 


.0 


.0 


.71 


-5 


13 


361 


.0 


.0 


.71 


-10 


14 


431 


.14 


.0 


.70 


-22.5 


15 


381 


.0 


.0 


.71 


-12.5 


16 


371 


.17 


.0 


.70 


-22.5 


17 


306 


.27 


.0 


.68 


-32.5 


21 


276 


.68 


.0 


.52 


-102.5 


25 


292 


.65 


.0 


.54 


-90 


32.5 


141 


.61 


.0 


.56 


-90 


41 


94 


.37 


.0 


.66 


-55 


90 


45 


.37 


.0 


.66 


-70 


139 


27 


.45 


.0 


.63 


-70 


166 


26 


.40 


.0 


.65 


-75 


250 


18 


.37 


.0 


.66 


-80 


400 


() 


.08 


.0 


.08 


-80 
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FIG. 1. Experimental spectra from the Orsay experiment [8] of the 4.438 MeV 7-ray line (open triangles) from the proton 
induced spallation of ^^O at Ep — 20 MeV and results of a Monte-Carlo simulation of the experiment (full line) with a constant 
differential cross section in the center-of-mass proton diffusion angle and a mean excitation energy of Ex = 14.75 MeV. Note that 
the experimental spectra at the different detector angles have been corrected for the individual relative detector efficiencies. The 
fact that the simulation with isotropic angular distributions predicts correct relative line intensities (with a common absolute 
normalization factor) indicates therefore an isotropic angular distribution of the 7-rays in the laboratory. 
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FIG. 2. (Upper part): Line shapes for the 4.438 MeV 7-ray at six different detection angles from proton inelastic scattering 
off ^^C at = 19.75 MeV. Triangles: count spectra of the Orsay experiment [8], full line: calculated line shapes with the 
population amplitude method. 

(Lower part): Laboratory 7-ray angular distribution from that reaction. Full squares: data from the Orsay experiment, full 
triangles: experimental data of the Washington experiment [11] at = 20 MeV, normalized to the Orsay data, full line: 
calculated distribution with the population amplitude method. 



14 




FIG. 3. Same as Fig. y, except full lines: line shapes and 7-ray angular distribution resulting from OM calculations. 



15 




20 40 60 80 100 120 140 160 180 



FIG. 4. Experimental differential inelastic cross sections ^'^C{a,a'y^C (filled squares, Refs. [31-37]) and results of OM fits 
(full lines). 
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FIG. 5. (Upper part): Line shapes for the 4.438 MeV 7-ray from proton inelastic scattering off ^^C summed for Ep — 
8.6-20 MeV. The summation has been done with a weight factor for each proton energy following the thick-target interaction 
probability of a incident proton spectrum as given in Eq. (p^. Parameters of the spectrum were s=3.3; -Eo=30. The summed 
experimental spectra from the Orsay experiment [8] are presented by open triangles while the summed spectra constructed 
with the population amplitudes method are shown by the full line. 

(Lower part): Summed laboratory 7-ray angular distributions. The summed experimental distributions are constructed from 
Legendre-polynom fits of the Orsay data and are shown by the full line, while the dashed line shows the angular distribution 
resulting from the population amplitudes method. 
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FIG. 6. Emitted line profiles of the 4.438 MeV 7-ray for flares with isotropic energetic particle distribution. Curves labelled 
1 and 2 are the profiles resulting from energetic proton and a-particle interactions with '^^C and '^^O with spectral parameters 
s—2.4; _Eo=300 MeV and s=3.3; Eo=30 MeV, respectively. Curves labelled 3-6 show the different contributions for the soft 
particle spectrum. They are in detail: (3) ^^C(p,p7), (4) ^'^0(p,pa7), (5) ^^C(a,Q7), (6) ^^0(a,2a7). 
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FIG. 7. Emitted line profiles of the 4.438 MeV 7-ray for downward flares at the solar limb, i.e. for a detection angle of 
90° with respect to the energetic particle direction. The upper two curves are for an energetic particle spectrum with s=2.4; 
£0=300 MeV, the lower two curves are for s=3.3; Eo=30 MeV. Full lines are the profiles calculated with population amplitudes 
from this work; for comparison, the dashed lines are calculated with population amplitudes from the parameterization proposed 
by Murphy, Kozlovsky and Ramaty [5] 
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FIG. 8. Emitted line profiles of the 4.438 MeV 7-ray for two fan-beam flares and an isotropic flare with energetic particle 
spectrum parameters s=3.3; i?o=30 MeV. The full line is for a fan-beam flare occuring at the solar disk, with an angle O/jare 
between the line-of-sight and the flare symmetry axis of 135° and an angular distribution proportional to sin®(0p), &p being 
the angle between the particle direction and the flare symmetry axis. The dashed line shows the profile for the same flare 
occuring at the solar limb (0/ioT-e=9O°), and the dashed-dotted line the profile for an isotropic interaction probability. 
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